


o 

0^ 



Romal n. 978 -1993 

ROM2F/93/38 

Revised version (May 1994) 



Deep Inelastic Scattering 
in Improved Lattice QCD 



^ ■ I. The first moment of structure functions 

> 



Stefano Capitani ^ 



^ ' Giancarlo Rossi ^ 

> 

k> i ^ Dipartimento di Fisica, Universita degli Studi di Roma "La Sapienza" , 

^ '■ and INFN, Sezione di Roma, P.le Aldo Moro 2, 1-00185 Roma, Italy 

^ Dipartimento di Fisica, Universita degli Studi di Roma "Tor Vergata" , 
and INFN, Sezione di Roma 2, Via della Ricerca Scientifica, 1-00133 Roma, 
Italy 



e-mail: 

capitani@vaxrom.romal.infn.it 

rossig@vaxtov.roma2.infn.it 



ABSTRACT 



1 Introduction 

Deep Inelastic Scattering (DIS) experiments, which started at the end of the 
sixties, shed hght on the properties of strong interactions and on the structure 
of hadrons (for a comprehensive review see [Q). The study of hard processes 
led to the formulation of the parton model and to the hypothesis that quarks 
are the hadronic inner constituents p|, |^. Further theoretical developments 
led to the discovery of asymptotic freedom 0, ^ and eventually to the formu- 
lation of Quantum Chromo-Dynamics 0. Today DIS experiments are still 
a very important tool of investigation for short distance physics phenomena 
and remain the best test of QCD as a theory of strong interactions. 

This paper is the first of two, focused on the computation of the renor- 
malization constants and mixing coefficients of the operators of rank two 
and three that are related to the first two moments of the DIS structure 
functions via the Wilson operator expansion of the product of two weak or 
electromagnetic currents. The knowledge of the hadronic matrix elements of 
the Wilson operators is necessary for the theoretical evaluation of the mo- 
ments of the x-distribution of quarks and gluons inside the hadrons. Our aim 
is to calculate to 1-loop renormalization constants and mixing coefficients of 
the operators related to the first and second moment in the framework of 
the nearest neighbor improved lattice QCD [0, ||, ^ |10|. These constants 
are needed to renormalize the lattice operators and be able to extract phys- 
ical hadronic matrix elements from numbers obtained in Monte Carlo QCD 
simulations. 

We present in this paper the general set up for these calculations and the 
results for the operators related to the first moment, leaving the discussion 
and the presentation of the results for the operator related to the second 



moment to a forthcoming paper |11|. There are two main reasons for the 



subdivision of the whole material in this way. One is that new high-statistics 
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Monte Carlo data for the first moment of DIS structure functions, produced 
in simulations employing the improved lattice QCD action, will soon be avail- 
able. The second reason, not unrelated with the previous one, is that, since 
simulations concerning the second moment have been given lower priority, it 
has not been yet decided what precise lattice expression for the correspond- 
ing rank three operator will be finally taken. This is not immaterial from our 
point of view: the perturbative calculation of the renormalization constants 
of the rank three operators are barely within the reach of algebraic manip- 
ulation programs such as Schoonschip or Form, because of the fantastically 
large number of terms generated in the course of the calculation. 

Renormalization constants and mixing coefficients of the rank two oper- 
ators have already been calculated some time ago using the standard Wilson 
action with not completely consistent numerical results fl^, |13|, |1^, [1^, |16|, |1^ . 



In this work, beside checking previous calculations, we present the results 
obtained for these constants by using the nearest neighbor improved QCD 
action (also known as the "clover- leaf" action), according to the improvement 
program started by Symanzik in ref. [p!8| . 

The use of this action has been proven ||^, |^, |^, |1^ to remove from on-shell 
hadronic matrix elements all terms that in the continuum limit are effectively 
of order a, a being the lattice spacing. Recent Monte Carlo simulations 
[0, 1^ show indeed a substantial reduction of the systematic errors related 



to the finiteness of the lattice spacing, though at the price of a slightly higher 
theoretical and numerical computational effort. 

Because of the great number of diagrams and of the complexity of the 
algebraic expressions involved in the calculation even in the case of the rank 
two operators, it has been considered necessary in this work to check all 
the results obtained by hand by means of suitable algebraic manipulation 
programs. 



The computation has also taken into account internal ferniion loops, hop- 
ing that the increasing power of dedicated computers would soon allow for 
fully unquenched QCD simulations. 

We have found discrepancies at the non-improved level, that is at the level 
of the standard Wilson action, with some of the numbers already published 
in the literature [O, O, |T1|, |T3[. They have both numerical and analytic 



origin. On the other hand, we are in complete agreement with the results of 



refs. [0] and [|T^ concerning the energy-momentum tensor. 

This paper is organized as follows: in Sect. 2 we discuss the general 
setting of the physical problem, and we immediately give the results for the 
renormalization constants and mixing coefficients of the relevant operators 
in a "ready-to-use" way. In Sect. 3 we present a brief introduction to the 
subject of improved lattice perturbation theory. In Sect. 4 the properties 
of the operators and their choice are discussed, in Sect. 5 we spell out all 
the necessary renormalization conditions, in Sects. 6 and 7 some peculiar 
aspects of the analytic and numerical calculations are explained, in Sect. 8 
we present the detailed results of our computation with a comparison with 
previous results, and in Sect. 9 we present some conclusions and an outlook 
of future lines of investigations. Some of the more technical aspects of the 
calculations can be found in the appendices A, B, C and D. 



2 General setting and results 

2.1 Moments of structure functions 

The DIS cross section is dominated by the field behavior in the hght-hke re- 
gion. The expansion of the product of, for instance, two scalar field operators 



on the light cone has the general expression |21 



A(x)5(0) ~Ec^.(^')^"^---^""O(^'L(0), (1) 






where the O^J^.'!^ are symmetric traceless tensors with vanishing vacuum 
expectation value: < 1 0|,^'.*|j^ 1 >= 0. The singularity of the coefficient 
functions is, up to logarithms, given by 

cnA^') ~ {x')^ ^ (2) 

and is governed by the twist [^ , 

r = do{N,i) — N, (3) 

of the operator, that is, by the difference between its physical dimension and 
the spin. The lowest twist operators are thus the dominant ones. 

In QCD the lowest twist operators appearing in the light-cone expansion 



of the product of two hadronic weak currents, relevant to DIS, are |2^, ^ 



25, m ra 



^Siv = T^^7[Mi H^, ■ ■ ■ D,^] (1 ± 75)17^ (4) 



Qgs = y Tr 

P 



^[fJ.\p -t^^2 ■ ■ ■ -t^MJV-l ^ f 



MJV-1 -^ Pf^N\ 



where the A-^'s are flavor matrices. The operators (^ are gauge invariant and 
all have r = 2; from now on we implicitly assume that they are symmetrized 



with respect to all Lorentz indices. S and NS superscripts refer to Singlet 
and Non Singlet flavor structures. 

In the unpolarized cross section the 75 contributions average to zero. The 
other contributions have matrix elements of the form 

< P\Ol!^^-,.N \P >= ^Nif^)p,,i ■■■Pm+ trace terms, (5) 

where /i is the subtraction point. They contain long distance non-perturbative 
physics, and are related to the general expression of the moments of structure 



functions by the equations p7 



|cixx^-i^fc(g2,x) = CN{qVfi')A^{li), (6) 

where we have put 

•^2 =f (7) 

with Fi, F2 and F3 the usual DIS structure functions as defined in [|I|. The 
Fourier transforms, Cn, of the Wilson coefficient functions are determined by 
the short distance perturbative behavior of QCD and are formally given by 

CN{q') = i<lTi^f Jd'x e'^^ c^ix'). (8) 

When renormalization effects are taken into account, a dependence upon the 
subtraction point /i appears in both sides of eq. (^). 

The Cat's are perturbatively calculable with renormalization group tech- 
niques. In case of mixing among different operators the Wilson coefficients 
satisfy a system of Callan-Symanzik equations 

^"1- + f^i9)4-)cN,i = T.i^N{9))ikCN,k, (9) 



dfi dgj 



k 



where the anomalous dimensions, '^n, have been arranged in a matrix 

{lN)lk = l^J- \og{ZN)Tk. (10) 

Renormahzation constants and mixing coefficients, {Z]^)ik, are defined by 
suitable sets of renormahzation conditions, allowing the construction of finite 
renormalized operators: 

ow)(/.) = E(^^(/^))/;cO(^''=). (11) 

k 

Equations (^ can be decoupled by diagonalizing the anomalous dimension 
matrix, so to have a new basis of operators that are multiplicatively renor- 
malizable. 



2.2 A summary of the results 

We have carried out, in the chiral limit, the 1-loop perturbative lattice cal- 
culation of the renormahzation constants and mixing coefficients of the fol- 
lowing rank two operators: 

0% = \^1[^.D,]^, (12) 

which is related to the first moment of the x-distribution of quarks, and 

0^, = ETr[F,,F,.], (13) 

p 

which is related to the first moment of the x-distribution of gluonsQ. The 

lattice expression of F^^ is given by 



^We are presently computing the renormalization constant of the operator 

Ol.^r^-^^llf.D.Dr]^', (14) 

which is connected to the second moment of the x-distribution of quarks pj| . 
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where Un,fiu is the usual plaquette 



'Jn,iJ,i' — '-^n,ij,'-^n+fi,u'-^n+u,fi^n,iy 



(16) 



The operators (|I2[) and (p!3[) are the first ones in the hst given above in 
eqs. (I). 

In the flavor Singlet sector, given the operators (0) and (0), the finite 
operators 0^^^ and 0^^^, with well defined anomalous dimensions, that are 
renormalized to their tree level values at p^ = yU^ [fia <^ 1), are given to 
1-loop in the full (unquenched) theory by the formulae (see Sects. 5 and 8 
for details) 



oU 



and 



^-Jz{N'^B,, + Nj{Bf^-B,, 



O^ 

^u 



^ + j-CF{Bgq-Bggj 



0% 
\17) 



6^ 

fMl/ 



Ctg ( 1 

— (72 log /ia + CpBqq + -NfBgg 






+ 



AC. 



as 



Air 



1 - -^ ( 72 log fia + ACpBgg + N.Bgg + NfBl 






(18) 



where N^ and Nf are respectively the number of colors and the number of 
flavors, and 



72 = 

Cp -- 


= {Nl-l)/2N, 


as = 


- 9l/^^- 



(19) 



The coefficients B for the case of the standard Wilson action and of the 
nearest neighbor improved action are reported (for r = 1) in Table 2.1. They 
can be derived from the more detailed Tables presented in Sect. 8. In Table 



Wilson 


Improved 


qq 


-3.165 


B-„ 


-15.816 


Bqg 


0.019 


B'n 


-1.041 


K 


-5.817 


B-„ 


-4.044 


(B^r 


-2.168 


(B'„y 


-19.425 


Bgg = -15.585 



Table 2.1 - Values of the constants B on the lattice for r = 1 



Bqq 


sails 

vertex 

self-energy 


-4 

8/9 

-1 


total 


-37/9 


Bqg 


-7/9 


Bgq 


-23/9 


Bgg 


11/12 


K 


- 23/18 



Table 2.2 - Values of the constants B in MS regularization 



2.2 we report for completeness the values of the constants B in the continuum, 
computed in the MS renormalization scheme [^6[ ^ . 

In the quenched approximation one has to put Nf = in the above 
formulae. Numerically, in terms oi (3 = 2Nc/gQ, for Nc = 3 one gets 



in the flavor Singlet sector 
i) for the Wilson action: 
9 



Oi. 



^7r2/5^^^ 






j^{K-B, 



27r2/3 



^^ fjLV 



(20) 



1 + 



1.776 



/? 






0.134 



P 






and 



fjLV 



1- 



37r2/3 
0.270 



log/za 



27r2/3 ^^ 






0.160 



/5 ^°^^"+ /3 
ii) for the nearest neighbor improved action: 

9 






(O 



1 \IMPR 
fiu) 



1 



17,^(3^'' 



O^ 

fif 



^ + ^j^i^i'i - ^i'^ 



1 + 



1.776 



/3 



O^ - 



1 + 



27r2/? 
0.596 



^Ol^MPK 



and 



(O 



2 \IMPR 



1- 



37r2/5 
0.270 



log fia — 



27r2/5 "^ 



(21) 



(O 



(J \IMPR 



(22) 



(O^J 



g \1MPR 



(23) 



log yua + 



0.801 



rQq ^IMPR 



fiu/ 



(3 -^ • (3 

where the explicit expression of (OJJ^)'^*^^^ can be found in eq. ( [50| ) of 
Sect. 6. 



9 



For the flavor Non Singlet operator 

, 1_ ^ A^ 

there is no mixing between O'^ and O^, and one gets 

1 






37r^f3 



log /ia 



27^^(3-^'' 



Qfq 



(24) 



(25) 



with obvious adjustments in notations in going from the Wilson to the 
improved case and with Bqq given in Table 2.1. 

3 Improved lattice QCD 

Lattice QCD represents today the only viable way of evaluating from first 
principles the matrix elements (^, needed for the computation of A]^ and, 
hence, of the moments of the structure functions. 

In this section we wish to briefly describe the improvement program of 
lattice QCD, as discussed in refs. [|T8|, ^ 0, |, g 0. 



The QCD Wilson action for one flavor / is, on a euclidean lattice [EU 



S^ 



LATT — (^ / . 



2a 



■^n+fj, 



+'^n+f,{r + l^^)U;t^^■^n\ + ^„ ( m/ + -- ) ^„ 



4r 



2 Z^ 



Tr 



'^ n,ii'^ n+n.u'^ ^ 



yfi'-^ n+ij,,v'-^ n+v,fi'-^ n.v 



f/J 



-iV. 



The gauge field A^^^ is introduced through the definition 

ig^at^Ai^^ 



Un„ 



(A 



1 A^^ 



(26) 



(27) 



Fermion fields are site variables, while the Un^fi^s are link variables. Un,fj, 
is associated to the link connecting the site n with the site n + fi. From 
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Wn+u 



^ n+u,fi 



^. 



n+fi+u 






^n 




U, 



n,ii 



^ri+uM 



^. 



n+fj, 



Fig. 1 - The ^z^-plaquette in lattice QCD 

this action one can derive the Feynman rules for lattice perturbation theory. 
Conventionally A„^^ will be taken to sit at the point n+fi/2, the middle point 
of the link {n,n + /i) (see Fig. 1). This choice greatly simplifies Feynman 
rules. 

The computations of physical quantities with Monte Carlo simulations are 
affected by statistical and systematic errors. The one which we are mainly 
concerned with here is the systematic error due to the finiteness of the lattice 
spacing. In fact for a renormalized finite lattice operator, Ol, we have the 
formal expansion 



P 



Or 



P 



Monte Carlo 



P 



o 



^■w.+^(") 



(28) 



where (p 



O 



p ) is the physical matrix element we want to extract from 

phys. 



Monte Carlo data and d is its physical dimension. We see from equation (Pq) 



that the possibility of estimating (p 



O 



p ) from Monte Carlo simulations 

phys. 



depends on the magnitude of 0{a) terms, which is generally between 20 and 
30 percent. The effect of improvement is the removal of all corrections that 
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in the continuum limit {qq ~ 1/loga) are effectively of order a, that is to 
say, of all terms that in perturbation theory are of the type 



a((/oT(loga)"~"a". (29) 

This result is achieved through the addition of "irrelevant" interaction terms 
to the Wilson action and by adjusting ("improving") the lattice expression 
of fermion operators, so as to cancel unwanted contributions in their matrix 



elements |18[, thus making faster the recovery of the continuum properties 
in the limit a ^ 0. It is estimated that with this method one can achieve a 
reduction of the systematic error due to the finiteness of the lattice spacing 



down to 5 - 10 percent [|1^, |2^ . 

The first proposal of improvement of the fermionic part of the QCD ac- 
tion presented, however, the drawback of having next-to-nearest neighbor 
interactions [^, |31], ^]. This feature is a problem in Monte Carlo simula- 



tions, given the calculational effort required for the numerical inversion of 
the fermion propagator (to date the only simulation carried out with the 
next-to-nearest neighbor improved action is that of ref. |3^). A great step 
forward was made by Liischer and Weisz, who introduced the notion of on- 



shell improvement [|28|: on-shell improvement is a weaker requirement than 
full improvement and easier to implement. Immediately after, Sheikholeslami 
and Wohlert proposed the nearest neighbor fermion action employed to- 
day and started the study of its properties. In ref. |^ the relation between the 
Sheikholeslami-Wohlert action and the next-to-nearest action was elucidated 
and it was shown that the Sheikholeslami-Wohlert action leads to on-shell 
improvement, with the consequence that all spectral quantities are free from 
terms that in the continuum limit are effectively of order "a" (eq. (^Qf)), 
provided that hadronic operators are at the same time appropriately "im- 
proved" . 
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The developments of these researches on on-shell improvement led to a 
practical recipe for its application to QCD [H, H, |10|. This consists in modi- 



fying the standard Wilson action by adding the nearest-neighbor interaction 
term 

^S{ = -igoO^ Yl -T '^n(^f^i^^n,fiu1pn- (30) 

n,fj,u 

Here Fn^^u is not the naive lattice "plaquette" 

1 
2igoa 



i^^n,fiv — ^._^ ^2 V "'Mi^ ^n,fj,u) {'^^) 



but rather the average of the four plaquettes lying in the plane /xz/, stemming 
from the point n: 






Because of the form of Fn^^u, the Sheikholeslami-Wohlert action is often re- 
ferred to as the "clover-leaf" action. 

The term (0), because of the antisymmetry of cr^j, matrices, does not 
alter the quark-gluon interactions with an even number of gluons. The in- 
teractions with an odd number of gluons are on the contrary modified. For 
1-loop calculations this results in the appearance of the new vertex 

iVy;ik, k') = ~go'-{t%. cos '^^tzJfk ^ ^^, sin a{k - k')^, (33) 

A 

where k and k' are the momenta of the incoming and of the outgoing fermions, 
and p is the Lorentz index carried by the gluon, to be added to the standard 
Wilson vertex 



{Vf;{k,k') = -g,{t\, 



a(k + k')p aik + k')p 
r sm h i-^p cos 



13 



(34) 



Besides this modification of tlie action, in tlie calculation of an n-point 
fermion Green function, one has to perform on each fermion field the rotation 
(see appendix A for the definition of the lattice covariant derivatives) 

ar 

y 

ar 



^ 



1 



^ 



v> 



(35) 



As we said, it was proven in ref. [^ that this recipe allows to get rid in on- 
shell matrix elements of all terms that in perturbation theory are of the form 



(pOD . The removal of these leading log terms leaves us with next-to-leading 
terms of the kind 



a(^o')"log' 



n-l 



(36) 



P 



O 



^W. + ^^^Z^^^"^ 



(37) 



which are effectively of order "a/ log a" . The improvement recipe thus lowers 
the difference between continuum and lattice from (^) to 

/ - \ IMPR. 

(p Ol p') = a 

\ I Monte Carlo 

Getting rid of also the next-to-leading terms would require a bunch of 
1-loop computations. They are necessary to adjust to next order the coeffi- 
cient in front of eq. (^0|) and the form of the transformation (|35|) . A crucial 
step forward in this direction has been carried out in ref. @] where the 
Qq correction to the factor in front of eq. (|30|) has been computed. This is 
enough to improve to next-to-leading log's the values of the hadron masses 
measured in Monte Carlo simulations p5|, BB . 



4 The operators 

We have computed, in the chiral limit, the four forward matrix elements 

<<l\OUq> <g,(r\Of^Jg,a> 



flV\ 



<Q\Oiu\q> <g,a\03\g,a> 



(38) 



IJ.V\ 
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where |g > is a one-quark state of momentum p and vanishing (renormahzed) 
mass and \g , a > is a one-gluon state of momentum p and polarization a. 

It should be noticed that beyond tree level there will be in general a mix- 
ing in the flavor Singlet sector between the quark operator O'^ and the gluon 
operator O^. However, in the "quenched" approximation {Nj = number of 
flavors = 0) one has < g,a\Ol^j^\g,a >= 0. In this case the matrix (^) 
becomes triangular and the operator O^ will not mix anymore with O^ (see 
eqs. (0), (|8D and (||) below). 

To write the operators (|12]) and (|T3D in a euclidean lattice, we need to 
take into account the changes in the symmetry properties of the Lagrangian 
[ P7] , |5B| due to the the breaking of the euclidean Lorentz group 0(4) down 
to the discrete hypercubic group 11(4): 

0(4) — > H{A) (39) 

(continuum —>■ cubic lattice). 

The hypercubic group H(4) is the group of the discrete rotations of the 
lattice onto itself. It is a finite subgroup of the ortogonal four dimensional 
group 0(4), consisting of 192 elements. H(4) admits 13 irreducible represen- 
tations. 

The breaking ( ^9]) implies that the operators belonging to irreducible 
representations of 0(4) may transform in a reducible way under H(4), thus 
inducing a mixing under renormalization among operators belonging to dif- 
ferent irreducible 0(4) representations. This mixing adds to the mixing be- 
tween (|12]) and (0), and causes further complications in the computation, 
even if the new operators have the same dimensions as the original one. If 
they have lower dimensions the coefficient of the mixing will inevitably con- 
tain negative powers of a, giving rise to power divergences: they can only be 
subtracted from the original operator in a non-perturbative way |^ . 
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To avoid all these unwanted mixings we need to choose particular values 



of the /ii/-indices in eqs. (|12D and (|13D , together with an appropriate lattice 
expression for the gauge field strength, F^^,, in eq. (|l^). 

First, we note that actually the quark operator has the explicit form 

0^, = iv^7[M/5.]^-5^.-tr, (40) 

and the like for the gluon case (0). The subtraction of trace terms, symbol- 
ically indicated by "tr" in eq. (0), reduces eq. @ to the simpler form 

< P\OiV-,^JP >= ^NP^, ■ ■ -P^^. (41) 

The operator O^j, with no trace subtraction is the sum of two irreducible 
representations of H(4): a four-dimensional one, corresponding to fi = u, 
and a six-dimensional one, corresponding to /i 7^ z/. With the subtraction 
of the trace the singlet is removed from the four-dimensional representation, 
which therefore becomes three-dimensional and cannot mix any more with 
the dangerous lower dimensional singlet operator ip ip with a hnearly diver- 
gent coefficient. However, large errors would result when the subtraction of 
the trace contribution is numerically performed on Monte Carlo data. The 
other option is therefore to choose /i 7^ z/. As we said above, this leads to 
an irreducible H(4) representation, so there will be no mixing with other 
operators. With this choice it will be necessary to give the hadron a momen- 
tum different from zero in one of the two fi or u directions. This causes a 
somewhat stronger sensitivity of Monte Carlo data to the granularity of the 
lattice. 



In actual Monte Carlo simulations [0, |14| the choice ^i ^ v has always 
been made, as the magnitude of the systematic errors associated with it is 
smaller than the one coming from the numerical subtraction of the trace 
contribution. In this paper also we will work with ^ ^ v. This, by the way, 
induces remarkable simplifications in the calculations. 
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We further restricted ourselves to the chiral case, that is to the case in 
which the (renormahzed) quark mass is zero. This choice corresponds to the 
hmiting situation of vanishing pion mass, to which Monte Carlo data are 
always extrapolated. 

As for the gluon operator, the same lattice approximation of the gluon 
field strength that is chosen in ( |52D turns out to avoid mixing of O^ with 
undesired operators [^, ^. In fact the operator (0) transforms like a 24- 



dimensional reducible representation and the product Pfj_p ■ Pp^, with fi ^ i', 
transforms like a 19 2- dimensional reducible representation. In particular 
Pfip ■ Ppu niixes with the operator ip ip with a power divergent coefficient. On 
the contrary, with the definition (|52D, which maximizes the symmetry of the 
expression of the field strength on the lattice, -F„,^i/ will transform like the 
direct sum of two three-dimensional irreducible representations correspond- 
ing to the two combinations E + B and E — B. Also in this case we choose 
11 ^ u to avoid the need for trace subtraction. 

5 Renormalization conditions 

The renormalization conditions connect the bare lattice operators on the 
lattice to finite operators renormahzed at a scale /i: 

^\^l) = Zlu{^la)0\a). (42) 

The constants Zi^ are fixed in perturbation theory by the same renormaliza- 
tion conditions used in the continuum. As we have discussed in the previous 
section, in the fiavor Singlet case there is a mixing between the quark and 
gluon operators (|12|) and (|l^) that have the same conserved quantum num- 



bers. However, with the choice ^i ^ v and the definition (|32|) of the gauge field 
strength we avoid mixing with other operators of possibly lower dimensions, 
and therefore the need for subtraction of power divergences. 
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We thus write: 



O" = Z,,0« + Z,gO^ 



O^ 



^aa^ + ^gg^ ; 



^gg^ 



(43) 



where the Z's are determined by imposing the renormahzation conditions: 



<q\d^{fi)\q> 
<g,a\d%n)\g,cr> 

<g|0%)|g> 
<g,(r\d3{fi)\g,a> 



<g|0^(a)|g>|;^.lV 



<g,a\0'^{a)\g,a>\';^ 



tree 



p-=^- 



<g|O»|g>|*^4>=0 



(44) 



In the perturbative computation of the two members of eqs. P^ ) we have 
in our hands the choice of the polarization index a of the external gluon 
(from which the Z's obviously do not depend). To simplify our successive 
calculations we decided to take cr 7^ /i, z/. 

In the tree approximation the amputated non vanishing matrix elements 
of the operators OJ^^ and O^j, are given (for /x 7^ z/, /x 7^ a, z/ 7^ cr) by 



<q\oiM\q>\X 



<9,(^\0Ua)\g,a> 



\tree 
lamp 



Ipt^Pu + (/i ^ Z/) 
-PliPy 



(45) 



From eqs. (|^) and ( pf ) the renormahzation constants and the mixing co- 
efficients of the rank two operators can be written for /ia <^ 1 in the form 

EH 
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Zqqifia) 

Zqg{^a) 
Zgqijja) 
Zggi^a) 



0:5/16 \ 

1 - —Cf ( y log fia + Bqqj 

~4^ (^y log/^« + ^OT 



47r 



Nf(-log^a + Bn+N,B^ 



'99 



(46) 



where Cp and as have been defined in eqs. (|T9|). 

The coefficients of the logarithms are the so called anomalous dimensions. 
They are physical quantities and obviously their values are not affected by 
the improvement procedure. The difference between the use of improved and 
non- improved action lies only in the finite constants B. Our results for the 
5's are reported in Sect. 8. 

If we want operators with well defined anomalous dimensions, i.e. opera- 
tors that are multiplicatively renormalizable, we must diagonalize the matrix 
( p!OD of anomalous dimensions. Its determinant, as it can be seen from 
is zero. The two eigenvalues are 



71 = 

72 = —CF + -Nf, 



(47) 



and the corresponding eigenvectors 



oi. 



and 



as_ 
47r 



{N,B,, + Nf{Bf^-B,,] 






+ 



0% 


= 


] 


Nf 


1 _ 


as 


ACf 


X 


4 


n 



«5^ 

47r 



0:5 / i 

^- ^-['l2^0g^ia + CpBqq + -NfBgg 

47r ^ 4 



(72 log fia + iCpBgg + N.Bgg + N f B ^^ 



[Bgq- 


Bqq) 


0'' 






(48) 


0^ 






'\ 


^^ IJ,U 




(49) 



The vanishing of one eigenvalue (71) of the anomalous dimension matrix 
means that there exists a conserved operator (O^). Physically O^ corresponds 
to the full energy-momentum tensor of the system. Notice also that in the 
quenched approximation {Nf = 0) the operator O^ is simply proportional to 
O'^, as no mixing is possible with O^ in absence of quark loops. 
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6 Details of the computation 

The improvement prescription discussed in Sect. 3 requires that in the fermion 
operators the spinor fields have to be rotated according to eqs. (|35D. This 
means that the improved operators are much more comphcated than the non- 
improved ones. For instance the rank two quark operator exphcitly becomes 



fQq ^IMPB^ 



fiuJ 



ar 






a'^r'^ 



16 



^ 



(50) 



(^^)7^D.^ ^-(V;^D,)7^ 

It should be noted that, although we are dealing with order a improve- 
ment, in calculating the Z's of eqs. (^3[) we have to take into account also 
the contributions formally of order a^ that arise when the rotations on ip e 
ip are both performed at the same time [^. These terms in fact do con- 
tribute to the constants B because 1/a^ divergences present in 1-loop lattice 
diagrams can compensate the factor a^ in front of them, thus giving rise to 
finite contributions to the Z's. Since in actual Monte Carlo simulations each 
ifj or ip field is subjected to the whole rotation ( pS] ) [ p!9[] , the full expression 
( |5nD must be considered in the calculation of renormalization constants and 
mixing coefficients. 

In the presentation of the computations in Sect. 8 we have for conve- 
nience separated all these various contributions. We will thus consider in 
turn diagrams with the insertion of the non-rotated (non-improved) opera- 
tor, diagrams with the insertion of only one rotation, either on tp or on ip, 
and diagrams with the insertion of both rotations. The operator ( ^0]) has 
been expanded in powers of ^^o up to order g^. This is sufficient for our 1-loop 
calculations. The corresponding formulae have been computed by hand and 
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then checked by means of the algebraic manipulation language Schoonschip, 
and are summarized in appendix A. 

Consequently to the use of improvement, there is a very large number 
of diagrams to be computed (see appendix B) and, a part from few cases, 
the algebraic manipulations of every diagram give rise to a huge number 
of terms. For this reason it has been necessary to check independently all 
the computations by using Schoonschip. We have to this end developed a 
general program that is able to automatically carry out all the algebraic ma- 
nipulations starting from the elementary building blocks of the calculation, 
represented by the expressions of propagators and vertices. The final output 
of the program is an analytic expression of the Z's under the form of a 1-loop 
integral. 

A large part of this work has thus been spent in developing an efficient 
Schoonschip code adapted to the problem at hand. The main difficulty re- 
sides in the fact that, although there are many built in instructions to deal 
with gamma matrices, Schoonschip has been conceived having in mind a 
continuum theory, which is invariant with respect to the (euclidean) Lorentz 
group. On the lattice, on the contrary, the theory is only invariant with re- 
spect to the hypercubic group, and unfortunately many simple and common 
terms like J2x l\Px sin ^a that arise on the lattice are not properly handled. 
Actually, the above term is wrongly reduced by Schoonschip to ]/sm.k\, be- 
cause two equal indices are by default assumed to be contracted. Furthermore 
the result depends on the order in which the various factors are encountered, 
because the Lorentz invariant summation is immediately performed between 
the first two indexed quantities encountered in the analytic expression. 

Being impossible to directly use the "gammatric" of Schoonschip as it 
is, a special routine has been developed to correctly treat gamma matrices 
on the lattice while using as much as possible of the built in Schoonschip 
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commands^. 

As we said, for each contribution the program starts with the expression 
of the corresponding Feynman diagram in terms of propagators and vertices. 
Before carrying out the Dirac algebra, one can simphfy the expressions by 
noticing that in the final loop integration one will have to integrate products 
of sines with functions of cosines and of H(4) invariant combinations of sines. 
Odd products of sines integrate to zero, while even products will lead to 
integrals of the type 

„ 2n 

X(/xi,...,/i2n) = d'^kf{cosk,^sin^ kx)Y[smki,^, (51) 

•' A 1=1 

where / is an H(4) covariant function. Exploiting the symmetry of the inte- 
gration measure under H(4), the integrand (|5T|) can be expressed in terms of 
a certain number of simpler integrals. For instance for n = 2 one has 

J(/ii, /i2, /Us, fJ'A) = / d^k /(cos k, Y^ sin^ kx) sin /c^, sin kf,^ sin fc^g sin fc^^ = 

/ d^k /(cos k, J2 sin^ kx) sin^ k^^ sin^ k^,.^ 6^^^^ 5^ 
•' X 

+ / d^k /(cos k, Y^ sin^ kx) sin^ k^^ sin^ /c^^ (5^^^3 5^ 
•' X 

+ / d*k /(cos k, Y sin^ ^a) sin^ A;^^ sin^ A;^^ 5^^^^ 5^ 
-' X 

- 2 ■ d'^k /(cos k, J2 sin^ ^a) sin"^ kf,^ ^mim2M3M4, (52) 

•^ A 

where 5^^^^^^^_^ is non-zero only if all the indices are equal. 

The case n = 1 is trivial. A 6-sine term {n = 3) instead gives rise to 



■'M3M4 



^M2M4 



^^2/^3 



a combination of 31 terms. In the calculation of < q\{01i^^^)^'^^^^\q > [[TT 



one also encounters some 8-sine terms (n = 4), each one giving rise to 379 
^We are willing to send our codes to anyone interested. 
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terms. A brief discussion of these expansions and of their derivation is given 
in appendix C. 

The rather comphcated Dirac algebra (we have products of up to 7 gamma 
matrices) is carried out at this point, exploiting the (5-functions that are 
expected to come out from the loop integration, as exemplified in eq. (|5^). 
The last step consists in extracting the Z factors by projecting the whole 
expression on the appropriate tensor structure (see eqs. (^5D). 



The CPU time needed to perform the whole analytic calculation varies 
according to the complexity of the diagram, going from 20 seconds up to 5 
minutes for the most complicated cases on a Sun 3 workstation. 

7 Loop integration 

The resulting analytic expressions must be finally integrated in momentum 
space over the first Brillouin zone. We want to compute all these integrals 
with a total error of about 1%, in order not to spoil the accuracy aimed at 
with the use of improvement. 

In almost all diagrams there will be an extra integration over a Feynman 
parameter a. Since it turned out to be exceedingly time consuming to eval- 
uate our five dimensional integrals with an error of less than 1%, we decided 
to perform analytically the a integration, writing the result in terms of the 
functions 

F„„(^/(^),^(A:)j = I da |y^^y^^^(^- (53) 

These functions, for the values of n and m relevant for this work, are tabu- 
lated in appendix D. The remaining four dimensional integrals are afterwards 
numerically computed with the most naive rectangular method, i.e. by sum- 
ming over a uniformly distributed net of points. 

All integrals are at most logarithmically divergent in the limit a/x — ^ 
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and have been computed by adding and subtracting to them known integrals 
with integrand functions having identical leading "1/A;^" and "1/A;"^" lattice 
behavior, as the original integrands. The final step consists in evaluating 
numerically the resulting finite integrals. In this way one can easily obtain 
the required precision with reasonable computing times. 

As for the "1/A;^" part, one subtracts and adds with the appropriate 
coefficients the integral 

'•^ d^k 1 



U(27r)4 [4EAsin^ 
whose value can be found in Table 7.1. 



Zq, 



(54) 



Zo 


0.15493339 


Zi 


0.10778131 


Fo 


4.36922523 


IE 


0.57721566 



Table 7.1- Values of some useful lattice constants 

For the logarithmically divergent part we actually use a double subtrac- 
tion. To explain the method, let us consider a typical integral of the form 



I = da 



N{k) 



7-^ (27r)4[L)(A;;a) + a(l-a)^V]2 



, N{0)=^0, (55) 



where 



D{k; a) = 4aJ2 sin' ^ + (^ " ") E ^in' kx + Ar^Y. ^in' ^)'] (^6) 



and, for instance, 



N{k) = cosk^. 



(57) 
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/ could be computed by summing and subtracting to it the known integral 

d'^k N{0) 



da 



-n (27r)4 [4 Ea sin^ f + a(l - a)fi^a 



,2^212 



N{0) 



167r2 



log /i^a^ — / da\oga(l 

JO 



a) 



+ 0{afj,), (58) 



7s + -^0 

where Fq is given in Table 7.1 together with the Eulero-Mascheroni constant 
'jE- In this way one can write 

iV(0) 



/=(/-/o) 



167r2 



1 2 2 

log/i a 



IE 



0{afi). (59) 



/ — Jo is finite as a/i —i> and can be safely evaluated in this limit. However, 



the denominators in (|55D and ([58|) are different, and this can lead to numerical 
errors larger than 1%, if one does not use a sufficiently large number of 
points, because the matching of divergent parts may not be sufficiently good 
in correspondence to the smallest value of k in the mesh of points one has 
taken. 

The results presented in literature for this kind of integrals generally use 
Monte Carlo integrations. With this method one can hardly explore the 
low k region with an adequate number of samplings, therefore there is a 
(small (?)) finite rest overlooked by the numerical integration. To overcome 
this problem (which is also present, though to a lesser extent, in the rectangle 
integration) and to reduce the number of different integrals to be computed 
numerically, we have employed a more refined subtraction procedure. 

The idea is to perform a double subtraction in order to match separately 
the form of the numerators and of the denominators appearing in the various 
integrals. To do this, we introduce the auxiliary quantity 

d^k N{{)) 



h 



da 



-TT {2'KY[D{k-a)Y 



(60) 
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and we write 

J = /o + (/i-/o) + (/-/i). (61) 

Of course both differences, Ji — Iq and I — Ii, are finite as a/x — i> and can 
be computed directly at afi = 0. Explicitly one gets 

lim /-/i = / da ,^ ,, ;'^ -V^ (62) 

a/^-o^ '' Jo i_^ (27r)4 [D{k;a)]^ ^ ' 

and 



lim (Ji - Jo) = iV(0) ■ fdaj 



[D{k;a)]^ [4EAsin 



2 fc;^]2 



[4 Ea sin^ ^]' - [4a Sa sin^ ^ + (1 - a) [Ea sin^ ^a + Ar\Ex sin^ ^)^]]^ 
[4aEAsin2 ^ + (1 - «)[EAsin^ fc^ + 4r2(^, sin^ f )2]]2[4Ea sin^ ^? 



We see that in eq. ( |62D the two integrands have the same denominator, while 
in eq. (^) the two integrands have the same numerator. The cancellation 
of the logarithmic divergences, between the two terms of eq. (|63D , has been 
made explicit in the last equality by the exact compensation of the leading 
l/fc^ behaviors of the two integrands. 

The double subtraction method we have explained reduces to only five 
the number of different types of integrals that are in the end necessary to 
express all the differences, similar to those appearing in eq. (|6lD, that arise 
in the diagrams one needs to compute. 

In the case of diagrams which contain only gluons, we have used the 
method suggested by Caracciolo, Menotti e Pehssetto [^ which allows to 
iteratively reduce all integrals to few basic integrals which are known with 
great precision. 
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8 Results 

The results of our calculations are summarized in the Tables 8.1 to 8.6 pre- 
sented in this section. The contributions coming from the standard Wilson 
action (that is, the non- improved results), those coming from the terms of 
order a in the improvement and those coming from the terms of order a? are 
separately shown. Contributions coming from different classes of diagrams, 
according to the classification given in appendix B, are also separately pre- 
sented. In Tables 8.7 and 8.8 we detail the comparison of our results with 
those of ref. [p!7| . 

In the quenched approximation [Nf = 0) Bgg and B^ disappear from 
eqs. (^) and (|49|) , and even in the flavor singlet sector the quark operator 
does not mix with the gluon operator, as the matrix (^) becomes triangular. 
In the full theory {Nf ^ 0, "unquenched" ) we have on the contrary complete 
mixing between the quark operator, O'^, and the gluon operator, O^. To 
compute in this case the renormalization constants and mixing coefficients 
we need also to evaluate the diagrams of Figs. 6 and 7 of appendix B, cor- 
responding to respectively the matrix element < g,o-\{0'^^,)^'^^^^\g,a > and 
the quark loop contribution to the matrix element < g, <y\0^^\g, a >. 

In Table 8.5 one can find the analytic expression of Bgg leading, through 
the use of Table 7.1, to the numerical values shown in Table 8.4. 

The contributions to < g, a\Of^Jg, a > of the diagrams labeled as "tadpo- 
le-QCD vertex" and "tadpole-quark loop" in appendix B are zero, and have 
not been inserted in Tables 8.4 and 8.5. 

As emphasized in the previous section, the accuracy of our results is bet- 
ter than 1%. The differences we have found by comparing our results for the 



Wilson case with results previously appeared in the literature [|12| |T3|, |1J, |15 
come from either insufficient accuracy in the Monte Carlo numerical integra- 
tion or trivial algebraic mistakes. However, we are in complete agreement 
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r 


Wilson 


0{a) impr. 


Ola^) impr. 


total 




0.2 


-1.970 


-1.746 


-0.026 


-3.742 




0.4 


-3.081 


-4.014 


-0.173 


-7.267 


SAILS 


0.6 


-3.935 


-5.524 


-0.487 


-9.946 




0.8 


-4.576 


-6.458 


-0.972 


-12.007 




1.0 


-5.077 


-7.041 


-1.626 


-13.744 




0.2 


0.555 


-0.132 


0.017 


0.440 




0.4 


1.038 


-0.470 


0.109 


0.678 


VERTEX 


0.6 


1.535 


-0.821 


0.290 


1.004 




0.8 


1.952 


-1.128 


0.555 


1.379 




1.0 


2.293 


-1.389 


0.904 


1.808 




0.2 


-6.102 


-0.611 





-6.713 




0.4 


-4.326 


-1.513 





-5.839 


i SELF-ENERGY 


0.6 


-2.762 


-2.311 





-5.073 




0.8 


-1.465 


-3.013 





-4.479 




1.0 


-0.381 


-3.646 





-4.027 




0.2 








0.006 


0.006 


operator TADPOLE 


0.4 








0.024 


0.024 


+ 


0.6 








0.053 


0.053 


i leg TADPOLE 


0.8 








0.095 


0.095 




1.0 








0.148 


0.148 




0.2 


-7.517 


-2.489 


-0.003 


-10.010 




0.4 


-6.369 


-5.996 


-0.040 


-12.404 


TOTAL 


0.6 


-5.161 


-8.657 


-0.143 


-13.961 




0.8 


-4.090 


-10.600 


-0.322 


-15.012 




1.0 


-3.165 


-12.076 


-0.575 


-15.816 



Table 8.1 - Values of 5, 



?<? 
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r 


Wilson 


0{a) impr. 


0{a^) impr. 


total 


0.2 


-5.184 


-0.029 


-0.010 


-5.223 


0.4 


-1.609 


-0.298 


-0.038 


-1.945 


0.6 


-0.704 


-0.582 


-0.063 


-1.349 


0.8 


-0.276 


-0.782 


-0.100 


-1.158 


1.0 


0.019 


-0.906 


-0.154 


-1.041 



Table 8.2 - Values of B, 



qg 



r 


Wilson 


0{a) impr. 


total 


0.2 


-3.722 


0.217 


-3.505 


0.4 


-4.372 


0.628 


-3.744 


0.6 


-4.938 


1.044 


-3.894 


0.8 


-5.413 


1.427 


-3.986 


1.0 


-5.817 


1.773 


-4.044 



Table 8.3 - Values of B, 



gq 





Wilson 


SAILS 


-4.453 


VERTEX 


5.019 


operator TADPOLE 


-40.024 +47rViV2 


i leg TADPOLE 


19.207 -2ti'^ /Nl 


i gluon LOOP 


2.274 


i GHOST 


0.198 


TOTAL 


-17.778 +27rVAr2 



Table 8.4 - Values of B, 



aa 
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SAILS 


192 97r2 + 24^'' 48^^ 24vr2 ' l^° ^^j 


VERTEX 


13 23 53 ^ 1 3 / \ 

192 + 487r2 144^°+3^^+167r2-r° ^V 


operator TADPOLE 


3 4 1 
64 3^°+4iV2 


1 leg TADPOLE 


32 12^'' 8Ar2 


1 gluon LOOP 


1 13 17 19 / \ 
32 ' 727r2 36^^ ' 1927r2 ' l^° ^V 


i GHOST 


11 1 / \ 


727r2 72^*^ ' 1927r2 [^' 'V 


TOTAL 


^ ^ ^ Z ^ Z + ^ 

64 487r2 16^'' 16 "' ' 8A^2 



Table 8.5 - Analytic expressions of Bgg/{lQiT'^). The constants Zq, Zi, Fq and 

^E can be found in Table 7.1. Each line in this Table equals the sum of the 

corresponding two in Tables 8.7 and 8.8. 
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r 


Wilson 


0{a) impr. 


total 




0.2 


-13.473 


-27.475 


-40.948 




0.4 


-5.643 


-32.141 


-37.784 


i quark LOOP 


0.6 


-3.512 


-27.057 


-30.570 




0.8 


-2.656 


-21.581 


-24.237 




1.0 


-2.168 


-17.257 


-19.425 



Table 8.6 - Values of B^g 



with the results for the energy-momentum tensor given in refs. |16] and [17 



as well as with the gluon self-energy reported in [JT]. 

To check the value of Zgg, it is enough to take the corresponding quantity, 
called Zi in Table 1 of ref. |]l6i, and multiply it by —8 and 167r^ to find 



Bgg (from ref. [16]) = -8 ■ 1671^ ■ (0.00252), 



(64) 



in agreement (within numerical accuracy) with the number given in the last 
line of Table 8.1: 



B, 



qq 



-3.165. 



(65) 



To compare the results for Zgg with those of ref. [r7| one must take into 
account that two different regularization procedures have been employed to 
deal with logarithmically divergent loop integrals. We have regularized all 
integrals by putting the external legs infinitesimally off-shell by an amount 
/i^ ^ l/a'^, while in ref. []1^ each term, J, is expanded up to second order 
in the four- vector of the external legs, p, and decomposed in H(4) tensors. 
Dimensional regularization is used to compute loop integrals. Calling J, as 



in ref. [17], the relevant terms of this expansion, it is easy to see that the 
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difference I— J (for small p) can in fact be computed as a continuum integral^, 
since its value only depends on the behavior of the integrand at fc = 0. 

To explicitly carry on the comparison with ref. Jl^ we have thus computed 



using dimensional regularization I — J for small p. We report their values 
in Table 8.7, apart from the logarithmic terms we always have singled out in 
all formulae. We collect in Table 8.8 the values of J, obtained from ref. JT7| 



remembering that in the Tables presented there our operator (|T^) corresponds 
in their notations to the tensor structure L3 — Liq. As announced, the sum 
J + {I — J) reproduce exactly the results given in our Table 8.5. 

Obviously physical quantities will be at the end independent from the 
chosen regularization procedure. In particular the dependence from the sub- 
traction point, /i, must disappear from physical hadronic matrix elements. 
In fact, consistently to each order in perturbation theory, the log fi terms get 
canceled in the product between the Wilson coefficients and the matrix ele- 
ments of the renormalized operators that are eigenvectors of the anomalous 
dimension matrix. The net result is that effectively the Wilson coefficients 
must be taken at a momentum scale a~^ and the operators (^8]) and (^) 



renormalized by reduced renormalization constants, obtained from the full 
expressions by dropping all logarithmic terms. 

The effective renormalization of the bare quark operator is rather small. 
Numerically for A^c = 3 at the typical values j3 = 2Nc/ Qq = 6 and r = 1, one 
gets in the quenched approximation for the unmixed quark operator 

Ol^ = 1.027 O;^^ Wilson case 

{0%y^^^^ = 1.134 {Ol^y^'^^ Improved case , 

where by the superscript ~ we mean that the log /i term contribution has 
been dropped from the expression of the Z's with the understanding that 
the corresponding Wilson coefficients are computed at a scale /i = 1/a. 
■^This is why only the genuinely lattice dependent quantity J is reported in ref. [ [L7[ . 
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SAILS 


7 /2 , , \ 
247r^ \e J 


VERTEX 


7 3 /2 , , \ 

0/1 2 + 1« 2 ■ +log47r 7£; 


operator TADPOLE 





i leg TADPOLE 





i gluon LOOP 


29 19 /2 , ^ \ 
144vr2 + 192vr^-U + ^°^'" ^') 


i GHOST 


1 1 /2 , , \ 


727r2 ' 1927r2 U ' °^ 'V 


TOTAL I- J 


73 

1447^2 



Table 8.7 - The contribution to the analytic expression of Egg/ilQir'^), cahed I — J 
in the text, according to the definition given in ref. [p^ ]. 
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SAILS 


1 7 31 19^ 7 /2^,\ 

^ H ^0 Zi ■ h Fo - bg 47r 

192 97r2 24 ^ 48 247r2 \ e ^ ) 


VERTEX 


13 3 53^ 1^ 3/2^, \ 
192 + 16.2 i44^°+3^^+16vr2-l i + ^° ^""^^V 


operator TADPOLE 


^ ^Z + ^ 
64 3^° + 4iV2 


i leg TADPOLE 


17^ 1 

32 12 '^^ 8iV2 


i gluon LOOP 


1 1 17^ 19 / 2 ^ , \ 
32 48vr2 36 ^° + 192.2 " ( " + ^o log4vrj 


1 GHOST 


-72^°+192.2-(-i + ^o-log4.) 


TOTAL J 


^ ^^ ^ Z ^7 + ^ 
64 18.2 16"'' 16 '^^ 8Ar2 



Table 8.8 - The contribution to the analytic expression of -63^/(16.2), called J 
in the text, according to the definition given in ref. [|1^]. 
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The gluonic contribution to the energy-momentum tensor operator, whose 
hadronic matrix elements are incidentally very difficult to measure in Monte 
Carlo simulations |]l3l, are somewhat larger. One gets in fact for O^ 



01^ = 1.296 09^ - 0.978 0«^ Wilson case 

(Ol^y^^P^ = 1.296 09^ - 1.099 {Op^^^^^ Improved case . 

There are suggestions that the large renormalizations in the gluon sector 
come from the choice of the bare lattice coupling as expansion parameter 
Actually, a very large renormalization factor relates the bare lattice 



coupling to the continuum one (as calculated for example in the MS scheme 
[0). It is also believed that the large contributions coming from tadpole 
diagrams are the main source of numerical discrepancies between continuum 
and lattice perturbation theory. 

According to ref. |4^, a resummation of tadpole contributions could 



be effectively achieved if the perturbative series is expanded in terms of a 
"boosted" coupling constant 

9 = -^: (66) 

where Uq{< 1) should be taken as the non-perturbative mean value of the link 
operator. A convenient measure of uq can be obtained from the plaquette 
expectation value (eq. ^l]) 

/I \l/4 

«o = ( o Trf/n,,/,,.) • (67) 

Here a ffist question arises whether or not tadpole diagrams should be re- 
traced from perturbative results, as they are effectively already taken into 
account in the redefinition (|66D. 

Another problem is that Uq could also be measured from the renormal- 
ization of the Wilson hopping parameter, by writing 
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The definition (§^), contrary to the one given by eq. (|67D, depends upon the 
choice of the ferniion part of the QCD actionQ. It should be observed here 
that, when using eq. (|68|), one should for consistency perhaps also introduce 
an explicit correction for the large deviation between the perturbative and the 
non-perturbative quark mass renormalization[|, besides changing the value of 
the effective expansion parameter. 

Because of the many subtleties involved in the whole question of "boost- 
ing" perturbation theory, we like to refrain from attempting any numerical 
evaluation of the effects of these corrections to our numbers, and we rather 
like to present the results of our calculations as they appears in terms of the 
bare coupling constant, leaving any other numerical consideration to actual 
applications. 

9 Conclusions and outlook 

As we have briefly discussed in Sect. 3, a simple modification of the fermion 
action, accompanied by a suitable redefinition of hadronic operators, leads 
to the theoretical expectation of an "effective 0(a)" improvement in on-shell 
hadronic matrix elements |p. This was indeed beautifully confirmed by the 
exploratory numerical Monte Carlo simulations of refs. |]T^ and 



After these encouraging results there has been a general agreement that 
the new high-statistics QCD Monte Carlo simulations that are planned to 
use next generation dedicated machines (e.g. APE with > 6 Gfiops) should 
all be performed with the "clover-leaf" improved action. 

This paper and the forthcoming one, that are dealing with respectively 



^It is interesting to note that, unlike the Wilson case, the two definitions (JCT) and (pi) 
seem to lead to numerically very close values in the case of the "clover-leaf" action Eaj . 



5 



We wish to thank G.Martinelli for a discussion about this point. 
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the first and the second moment of DIS structure functions, are part of the 
general project aimed at the perturbative recalculation with the "clover-leaf" 
improved action of all the renormalization constants and mixing coefficients 
of the fermion operators whose matrix elements are phenomenologically rel- 
evant and susceptible of being measured in Monte Carlo simulations. 

To complete the analytic part of the above program only the mixing 
coefficient between the dimension six four- fermion A/ = 1/2 effective weak 
Hamiltonian and the dimension five operator appearing in eq. (^) is missing. 
In the standard Wilson theory this finite (thanks to the GIM mechanism) 
coefficient was computed in ref. [^. Its calculation in the case of the "clover- 



leaf" improved action is presently under way 

We think it is fair to conclude this paper with a remark on the reliability 
of lattice perturbation theory itself, which looks even more problematic than 
continuum perturbation theory. For instance, gluonic corrections tend to be 
anomalously big because of large tadpole contributions. As we remarked in 
Sect. 8, one possible remedy could be to use a better expansion parameter 
than the bare coupling constant, as advocated in refs. |^ and ^^, in a way 



to effectively incorporate a resummation of these large tadpole diagrams. 
A perhaps more promising and consistent way to proceed has been, how- 



ever, recently suggested in ref. [^. The idea is to extract renormalization 
constants directly from Monte Carlo data, by measuring matrix elements of 
hadronic operators between quark states and exploiting information coming 
from Ward identities. Such a procedure would eliminate altogether lattice 
perturbation theory from the game. This fact might be regarded as rather 
satisfactory for a numerical approach like the lattice one which, taken in its 
most extreme formulation, should refrain from using any kind of (approxi- 
mate or unreliable) perturbative calculations. 
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Appendix A 

Perturbative expansion of the vertex operators 

(0« J»'™ and Ol, 

In this appendix we give the perturbative expansion of the vertex operators 
(O^y)'^*^^-^, given by eq. (|50|), and O^^, given by eq. (|13D with the insertion 
of the definition (p^). We used the covariant derivatives defined as 

^/. V'n = — \Un,v'4^n+^, - ^,|-^,^,^n-J (69) 



^ 2a 
and the conventions 



n+pi^ n,u r n— /i^J^— MiA* 



^-^^=j4^.''^'^''^^^"M'i) (70) 



d'^q 



5?„ = /j^ €-'""?(,), (72) 

where the integrals are performed over the first Brillouin zone. Throughout 
this appendix external and loop momenta are expressed in lattice units. 

Since the full Fourier transform of the operator {O'j'^^)^'^^^^ is very com- 
plicated, we give here only the form it effectively takes when inserted in the 
diagrams of Figs. 2, 3, 4 and 7. Calling p the external incoming and out- 
going momentum, and k the fermion loop momentum, one finds, separating 
the various contributions according to their naive order in a: 

a) tree level 
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+ 



ar 



ia^r"^ 



2 J-n (27r)4 



^'^ i^ikf^^'^^m 



(73) 



-TT {27iy 



- ,^ sin A;a sin /ca' sin fc;, 

^(^)z^7a7m7a' V^(«;)- 



A,A' 



b) order go 
{0%) 
agor 



\IMPR, 



n 



ISO 



d^k 






^9Q_ 
2 

cosf 



A; + p 



V'(p)7a« 



cos( )^A^(p-/c) ^(/c), 



^ sin^A sinpA\ . . ,. 

7m7a + 7a7m ]Au[p- k) 



+ COS 



ia'^gor^ 



2 
d'^k 



, / sin pi, sin fcj, 

)a(7.7a^ + 7a7.^ 

A; + p 



TT (27r) 



:V^(P) Zl 7a7m7a' 



+ cos I 



A,A' 

k + p . 



COS 



]Ax{p-k) 
)^A^{p-k) 



sinkx' sin Pa 



xAx{p-k) 



sin /ca' sin k^, 



+ cos ( ^7^) A' A' (P - k) 



a a 

sin Pa sin p^, 



Hk). 



(74) 



~ 2 ■ ■" a a 

This formula is given in the kinematical configuration in which the incoming 
gluon momentum lands on the incoming quark leg. If the gluon is attached 
to the outgoing quark leg, one must exchange p and k. Notice that the first 
lines in eqs. ( |75D and (^) correspond to the non-improved expression of the 
operator. 

We do not give here the expression of the 0((7q) terms because, besides 
being extremely complicated, they are not actually necessary for our compu- 
tation. In fact we need them only when either the gluon or the quark legs are 
contracted to make a tadpole loop, and in this situation the tadpole directly 
f actor izes out. 

On the gluon operator O^ there is no effect due to the improvement. 



because the transformation (|35D acts only on spinors, and one finds: 
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a) tree level 



0^.(^ = 0)1*. 



1 



2^] (27r)4 



^^1^^2^(4)(^^ + ^^) 



- i^liQi) cos -^ sin qip - A^(gi) cos -^ sin q^^ 
■ - \Ap{q2) cos — sing2i. - A^{q2) cos — sing2p 



(75) 



b) order qq 



r^q / nM '''do rabcSr^ f dqi dq2 dq^ u) . , 

0^.(^^ = 0)150 = -^/ Ly (2^X8 ^^^91+^2 + ^3) 



1 

a 



A)^(gi) cos -^ sin gi^ - v4^(gi) cos ^ sin gi^ 



Al{q2)A: 



' <i2p {q2 + 2q3)p^ 
cos — - — COS 

2 2 



^ 931/ {q3 + '2q2)u' 
COS COS 

2 2 



(ga + 2g3)p (gs + 2g2)i. 
-2 cos cos 



Ab/ \ AC/ \ ■ ('?2 + q3)p ■ 

- A {q2)A {qs) sm sm gg^. 



+ 4(g2)Ae(g3)sin. 



?2 + qsju 



■ sm gsp 



+ (/x ^ u). 



(76) 



As for the OIqq) terms we do not report here their expression for the 
same reasons explained after eq. (|73|). 

Appendix B 
Diagrams 

In what follows we show all the 1-loop diagrams that have been calculated 
in the context of this work. Many of them have a label that allows an easy 
connection with the Tables presented in Sect. 8. 
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Actually, each graph in this appendix corresponds to several diagrams in 
the perturbative expansion. For example, the vertex correction of Fig. 2 




Fig. 



The graph that symbolically represents the 1-loop correction to the 



insertion of the {O'^^) operator. The insertion is indicated by a dot. The 

wavy line is a gluon. 



symbolically represents the sum of 12 different diagrams. They are all shown 
in Figs. 3.1 to 3.12, where the improved quark-gluon vertex ( ^3]) is repre- 
sented by a dot and the insertion of 0{a) and Ola"^) corrections to the vertex 
operators (see also eqs. ( |5I]| ) and ([THD) are respectively indicated by an open 
circle and a square. 





3.1 



3.2 
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3.3 



3.4 





3.5 



3.6 





3.7 



3.8 
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3.9 



3.10 





3.11 



3.12 



Fig. 3 - The 12 diagrams that contribute to the graph of Fig. 2. We have indicated 
the improved quark- gluon vertex with a dot and the insertions of 0{a) and O^a^) 
corrections to the vertex operators with respectively an open circle and a square. 

With this understanding we show below only the general patterns of the 
various graphs, as representatives of the full set of improved diagrams. 

The graphs needed for the 1-loop computation of the matrix element 
< Q\{0'li^)^'^^^\q > and the graphs contributing in the quenched approxima- 
tion to < g,(7\0fi^j^\g,cr > are shown in Figs. 4.1 to 4.8 and Figs. 5.1 to 5.13 
respectively, where the insertion of the vertex operator is always indicated 
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by a dot. 





4.1 - Vertex 



4.2 - Sail 





4.3 - Sail 



4.4 - Operator tadpole 




4.5 - Self-energy 




4.6 - Self-energy 
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4.7 - Leg tadpole 



4.8 - Leg tadpole 



Fig. 4 - The different types of graphs contributing to the 1-loop approximation 



of the matrix element < q\{Of^ 



\IMPR\ 



fll/J 



>. 





5.1 - Vertex 



5.2 - Sail 





5.3 - Sail 



5.4 - Operator tadpole 
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5.5 - Tadpole (QCD vertex) 5.6 - Gluon loop 





5.7 - Gluon loop 



5.8 - Leg tadpole 





5.9 - Leg tadpole 



5.10 - Ghost 
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5.11 - Ghost 



5.12 - Ghost 




5.13- Ghost 

Fig. 5 - The different types of graphs contributing to the 1-loop approximation of 
the matrix element < 5, a\0'j^^j\g^ a >, in the quenched approximation. The dotted 
line represents a ghost loop. 





6.1 - Quark loop 



6.2 - Quark loop 
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6.3 - Tadpole- quark loop 



6.4 - Tadpole-quark loop 



Fig. 6 - The different types of graplis contributing to the 1-loop approximation of 
the matrix element < g ^ (j\0^^\g , a > to be added in the full (unquenched) theory 
to the diagrams of Figs. 5. 



In the full (unquenched) theory there will be mixing between the flavor 
Singlet operators {O'^^)^^'^^ and O^^. To compute the mixing coefficients 
we have 

1) to add to the diagrams of Figs. 5 those of Figs. 6, in which a quark 
loop is present. 

2) to evaluate the off-diagonal matrix elements < g,a\{0'^j^)^'^^^\g,a > 
and < q\0^ \q >. The corresponding diagrams are listed in Figs. 7 and 
Figs. 8 respectively. 
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7.1 



7.2 





7.3 



7.4 




7.5 

Fig. 7 - The different types of graphs contributing to the 1-loop approximation 
of the matrix element < g, o'l(0^jy)^*^^'^|5', o" >. 
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^.1 



8.2 



Fig. 8 - The different types of graphs contributing to the 1-loop approximation 
of the matrix element < q\03 \q >. 

Appendix C 
Products of sines 

At some intermediate stage of the computation there appear expressions 
whose general structure is of the kind 

. 2n 

X(/ii,...,/i2n) = d'^kf{cosk,'^sin^ kx)Y[smkf,^, (77) 

By exploiting the H(4) covariance properties of the integrand, the integral 
(|77D is first reduced to a sum of simpler terms that have the form 

/ d'^kf {cos k, J2 sin^ kx) J] sin^"" k^, (78) 

•^ A M=l 

where n^ are integers ranging between and n and satisfying 

4 



J2n^ = n. 



(79) 



M=i 
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■'M3M4 



Ml7^M3 



For example in the n = 2 case one gets 

/ d'^k /(cos k, Y^ sin^ kx) sin^ /c^^ sin^ k^,^ 6^^^^ 5 
•^ A 

+ / d^k /(cos /c, ^ sin^ /ca) sin^ /c^^ sin^ k^^ 6^,^^^ 6^ 

■' A 

+ / d^k f {cos k, Y^ sin^ kx) sin^ fc^, sin^ k^^ 6^^^^ 6^ 
•' A 

+ / d^/c /(cos A;, ^ sin^ /c^) sin'^ A;^^ 5^^ 



^t^2^J■4 



^M2M3 



Ml7^M2 



^17^/^2 



1M2M3M4' 



(80) 



where S^^^^^^^^ is non-zero only if all the indices are equal. 

In the form (^Dj) this equation is not suited for the further algebraic ma- 
nipulations that Schoonschip will have to perform. To overcome this problem 
we have simply to implement algebraically the various conditions /ij 7^ ^j in 
eq. (|80| ) by writing 



sin^ k^^ sin^ k^^ 



^J■l¥=^J■2 
and the like. This leads to the formula 



sm fc^j sm fc^2 sm k^-^ "/ii/12 



(81) 



T(/ii, /i2, /U3, /i4) = / d^kf {cos k, Yj sin^ kx) sin /c^^ sin fc^j sin fc^g sin A;^ 
/ d^A; /(cos k, Y sin^ A;a) sin^ A;^^ sin^ A;^3 (5^^^2 5 



■'A'3M4 



-F / d^A; /(cos k, Y sin^ A;a) sin^ k,,^ sin^ A;^2 (^^,^3 ^^^^^ 
■' X 

+ / d^A; /(cos A;, ^ sin^ A;a) sin^ k^^ sin^ A;^^ 5^^^^ S^^^.^ 
•' X 

— 2 ■ / (i^'A; /(cos A;, ^ sin^ A;a) sin'* A;^^ 5^^ 



M2/i3M4- 



(82) 
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If n is higher, the way to get useful expressions is the same, but it is 
necessary to use relations of the kind (^Tj) several times. The expansions for 
n up to 4 are summarized in Tables CI to C.3. 



Type of term 


Number of permutations 


Weight of each permutation 




3 


1 


"M1M2M3M4 


1 


-2 



Table C.l - Numerical coefficients for the expansion oi Z{fii, fj,2, fJ-3, fj.4) 



Type of term 


Number of permutations 


Weight of each permutation 


"A A A " 


15 


1 




15 


-2 


"Ml/'2/i3M4/i5M6 


1 


16 



Table C.2 - Numerical coefficients for the expansion of X(/ii,/X2,/X3,/i4,/U5,/X6) 



Type of term 


Number of permutations 


Weight of each permutation 


"A AAA" 


105 


1 


"M1M2"M3M4"M5A'6A'7A'S 


210 


-2 


"A«1M2"M3M4M5A'6A«7A'8 


28 


16 


"Ml/i2/i3M4"A'5M6/i7M8 


35 


4 


"/ilM2AI3/i4M5M6M7M8 


1 


-272 



Table C.3 - Numerical coefficients for the expansion of 

I{fll, /U2, /X3, /i4, /"5, /"6, /^7, /"s) 
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A check on the correctness of these coefficients may be obtained by con- 
sidering the case in which in eq. ( [77|) all indices are equal. In this case there 
is only one contribution. This means that in each Table, if we sum the num- 
bers obtained by multiplying the weight of each permutation by the number 
of permutations, we must get 1, as it is immediately checked in all cases. 

Appendix D 
a-integration 

To reduce the computing time in the numerical integration of the integrals 
that involve fermion propagators, we have chosen to perform analytically 
the integration over the Feynman parameter a, and to leave for the numer- 
ical integration only a four-dimensional expression. To this end we need to 
compute integrals of the form 

^w-^w) = i'"° i/w+°Ig(tr - (83) 



i^\ fx — A . o , . o /x — A . 2 A\2 



In our calculations we explicitly have 

g(k) = 4Esin^^-[EsinHA + 4r^(Esin^T 

A ^ A A ^ 

f{k) = \Y^sm'k, + Ar'{j:sm'^r]. (84) 

A A ^ 

The functions F„m satisfy the recurrence relations 

d 

TTF^ nm = ~f^ ' ^ n m+1 (85) 

d 

— F„„ = -m ■ F„+i m+i (86) 
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which make simpler their computation and allow for a check of the formulae 
given below. 

The formulae needed in this work are: 



Fo2 = 


f-{9 + f) 












Fl2 = 


1 
9 


-log 1 + ^ — 

9 V fj 9 + f. 




F22 = 


1 
7' 


1-2^ log 1 + ^ +— ^ 
9 V fj 9 + f. 




F32 = 


1 

7' 


2 + 3 log 1 + 
.2 ^ ^2 V // 


/^ 1 




^■(^ + /)J 




Fno - 


9 + V 




-T 03 — 


2P ■ {9 + f? 


Fl3 = 


1 




V-{9 + f? 


F23 = 


1 


[1, /^^ , g\ 3(7 + 2/ " 




F33 = 


1 


9 \ fj 2-(^ + /)2_ 




F43 = 


1 

7' 


"1 0/ , rf\ (. , ^\ .2 7(? + 6/ ] 

3 + 6 log 1 + / . , .X2 




F53 = 


1 


"1 3/ , J-^ ^J\ (^ , g\ , ^3 9(^ + 8/ ] 


Fo4 = 


9' + 3^/ + 3/2 
3/3 ■ {9 + ff 


Fi. - 


9 + ^f 




r 14 — 


Qp ■ {9 + fr 


Fn. - 


1 




-^24 — 


3/-(^ + /)3 


F34 = 


1 


-log 1 + ^ -^ 
9 \ f) 6 




^ + 6/2 
/)^ 









(87) 
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